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We investigate the effect that the intrinsic spin-orbit and the inter- and intra-layer Rashba inter- 
actions have on the energy spectrum of either an unbiased or a biased graphene bilayer. We find 
that under certain conditions, a Dirac cone is formed out of a parabolic band and that it is possible 
to create a "Mexican hat" -like energy dispersion in an unbiased bilayer. In addition, in the presence 
of only an intralayer Rashba interaction, the K (K') point splits into four distinct ones, contrarily to 
the case in single-layer graphene, where the splitting also takes place, but the low-energy dispersion 
at these points remains identical. 
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I. INTRODUCTION 

The influence of spin-orbit (SO) interactions on a sin- 
gle layer of graphene is well known. Kane and Mele^l 
were the first who showed that the intrinsic SO (ISO) in- 
teraction not only can open a gap, but it also gives rise to 
a quantum spin Hall phase, due to localized edge states. 
On the other hand, the extrinsic Rashba SO interaction 
acts in the opposite direction, and tends to close the gap. 
Later, it was found that the Rashba interaction has in- 
teresting effects on its own, leading to a splitting of the 
Dirac point into four identical points P This splitting is 
missed in the low energy calculations used in Refs. 11121 

The presence of a gap controlled by the ISO interac- 
tion seemed promising, but it turned out to be smaller 
than originally expected by Kane and MeleP^The value 
of its coupling constant is still controversial, but is ex- 
pected to be in the range 0.0011 - 0.05 meV, 6 7 thus 
very small. The Rashba coupling instead, can be tuned 
to much higher values. For typical values of an exter- 
nal electric field (50 V/300 nm), the Rashba coupling is 
less than 1 meVP^ The effect of impurities can increase 
this value to 7 meV.^ However, recent experiments on 
epitaxial graphene grown on a Ni(lll) substrate showed 
that the Rashba coupling can reach values up to 0.2 eVp^ 

Shortly after the discovery of graphene, it was ob- 
served that bilayer graphene also exhibits remarkable 
phenomena. In bilayers, the low energy excitations are 
no longer Dirac fermions like in graphene, but massive 
chiral fermions In addition, bilayer graphene turns 
out to be a semiconductor, with a gap that can be 
tuned via a chemic al doping^ or by an external gate 

voitagePnannaig 

Although there are several studies for the effect 
of SO interactions in graphene accounting for differ- 
ent boundaries 19 ? (zigzag and armchair) and electron- 
electron interactions no investigations of the SO ef- 
fects have so far been performed in bilayer graphene, to 
the best of our knowledge. In this paper, we incorpo- 
rate SO interactions in a bilayer graphene system in the 
presence and absence of a bias voltage. We distinguish 
between the ISO interaction, which respects the lattice 
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FIG. 1: Relabeling of the graphene lattice. The length of 
the lattice vectors shown equals a, while the lattice spacing is 
a/y/3. Figure extracted and modified from Ref. l25l 



symmetry and the extrinsic Rashba interaction, which 
is only present if the lattice symmetry is broken. Here, 
we break the lattice symmetry by introducing an electric 
field. Depending on the orientation of this electric field, 
there can be both intralayer as well as interlayer Rashba 
interactions. Furthermore, we consider how the energy 
spectrum deforms if the layer is biased with an electrical 
voltage. 

As a main result, we find that the intralayer Rashba 
coupling in an unbiased bilayer not only splits the K 
(K f ) points into four, in a different way than it does 
for monolayer graphene, 3 but also creates a Dirac cone 
out of a parabolic band. In addition, we show that a 
fully spin-polarized Mexican-hat band arises in the en- 
ergy spectrum of an unbiased layer, purely due to SO 
interactions. 

The Mexican-hat-like dispersion appears in a variety 
of physical systems. This kind of spectrum, with a line 
of degenerate low-energy points forming a ring, was first 
discussed by Brazovskii,^ who showed that it leads to a 
"weak" crystallization transition. In cold atoms physics, 
a Mexican-hat-like dispersion appears and gives rise to 
topologically different ground states in SO Bose-Einstein 
condensates.^ The Mexican hat is known in high energy 
physics as well, where, for example, the Higgs mechanism 
is expected to be responsible for the mass of the vector 
bosons. For bilayer graphene without SO interactions, 
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FIG. 2: (Color online) Behavior of a graphene sheet for dif- 
ferent values of the SO parameters, (a) Aso = 0, £r = 0. (b) 
A so = 0.5, t R = 0. (c) Aso = 0,t R = 0.2. (d) A so = 0.5, 
t R = 0.2. In all figures t = 1, a = 1, A; y = 27rv/3/(3a). 



the energy dispersion has the Mexican-hat shape if the 
layers are biased. In this case the dispersion was shown to 
potent ialize electron-electron interactions, thus leading 
to a ferromagnetic instability. 24 

The outline of this paper is the following: To render 
the comparison with a monolayer sample easier, we recall 
some results for monolayer graphene in S ec. \U\ Next, we 
set up a model for bilayer graphene in Sec. |III[ after which 
we include SO interactions within the layers in Sec. [IV| 
We also include Rashba interactions between the layers. 
This is done in Sec. [V] Finally, we draw our conclusions 
in Sec. ED 



II. SO INTERACTIONS IN MONOLAYER 
GRAPHENE 



In graphene, the carbon atoms arrange themselves in 
a honeycomb lattice. Because there are two inequivalent 
positions for the carbon atoms, this honeycomb lattice 
can be seen as a triangular lattice with two atoms per 
unit cell, called A and £?, see Fig. [I] In the tight-binding 
approach, one assumes that the electrons are localized 
around the lattice sites and that they can hop from one 
lattice site to the next (nearest neighbor hopping). The 
noninteracting Hamiltonian is then given by 



H = t J2 [4(H*)MRi + «j) + fc.c] , (1) 

J=l,2,3 



where a is the spin index, i runs over the A-sublattice 
sites and j over the nearest neighbor vectors, which with 
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FIG. 3: (Color online) Splitting of the Dirac cones due to 
the Rashba interaction. For A so = no gap opens. In this 
figure: t = 1, a = 1, k y = 27ry / 3/ (3a), Aso — 0, and t R — 0.2. 



a lattice orientation as in Fig. [T] are denned by 



^(0,1). 
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The constant t is the hopping parameter 2.8 eV) and 
a/V3 « 0.142 nm is the lattice spacing. 

After relabeling the lattice as in Fig. [TJ one can bring 
the Hamiltonian ([T]) into the form 



H = t 



aj.(m, n)6cr(m, n) + a£(ra, n)b a (m — 1, n) 



+ a£(ra, n)6cr(m, n — 1) + h.c. 

By performing a Fourier decomposition, the free Hamil- 
tonian reads 

Ho = tj d 2 k^(k)Ml xA ^(k), (2) 
where ^t(fe) = (at(fe) T ,at(fe) i ,6t(fe) T>6 t(fe) i ) > 



M 4 ° x4 = 



/ 7k \ 

7 k 

7k 
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and |7k| 2 = 3 + 2 cos(ak x ) + 4 cos(ak x /2) cos(V3/2ak y ). 
The eigenvalues of Hq are the well known energy bands 
of graphene,^ E± = ±t|7k|. Both bands are degenerate 
with respect to the spin degrees of freedom. 

If one wants to include SO interactions into the 
graphene system, one has to distinguish between the ISO 
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FIG. 4: (Color online) Energy spectrum of bilayer graphene for different values of the SO interactions. In this figure the layer 
is unbiased, hence V = 0. Other parameters have the values t = 1, t± = 0.2, a = 1, and /c y = 27r/(v / 3«). 



interaction and the extrinsic Rashba term. The ISO in- 
teraction does respect all symmetries of the graphene lat- 
tice and has the form of a next-nearest neighbor (n.n.n.) 
hopping term, 



Hso = itso ^ 



In this expression, tso is the n.n.n. hopping amplitude, 
s z is the z-component Pauli matrix describing the spin, 
and cj is either aj or 6j, depending wether the index j 
labels an A or B-sublattice site, respectively. The factor 
Uij is +1 if the n.n.n. hopping is anti-clockwise and —1 
if it is clockwise. Note that this term describes hopping 
within the same sublattice. 

Using the relabeling of the lattice shown in Fig. [I] and 
performing a Fourier decomposition, the ISO Hamilto- 
nian can be rewritten as 



where 



so 

4x4 



(j) 

-r] 

-rj 

\0 Or] 



(3) 



r] = [1/(3^3)] [2sin(ak x ) - 4sin (ak x /2) cos (y/3ak y /2)] 

and A S o = 3^3 t S o- 

The extrinsic SO interaction is the Rashba term, which 
is only present if the lattice symmetry is broken. This can 
happen if the graphene sheet couples to a substrate or if 
an electric field is present. For a perpendicular electric 



field, E = Ez, the Rashba coupling has the form of a 
nearest neighbor hopping term and is given b ^ 2 * 3 * 26 * 



Hr = it R ^ 4 ( S X ^ij) ' ZC 3 + h ' C 'i 



(4) 



<*,j> 



where the hopping amplitude tR can be tuned by chang- 
ing the electric field strength, s is the vector of Pauli 
matrices and d^- is the unit vector that connects the i 
and j lattice sites. This term describes nearest neighbor 
hopping, but it only couples nearest neighbors with op- 
posite spin. This is clearly seen if we rewrite this term 
in the same way that we rewrote the other terms, 



where we have defined 



t R / d 2 k^(k)M* x ^(k), 



(5) 
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The total Hamiltonian can be obtained by collecting Eqs. 
0, (§, and (§, 



H 
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In Fig. [2] we show the behavior of the energy disper- 
sion for a graphene sheet for different values of the SO 
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FIG. 5: (Color online) Due to the Rashba coupling, the spin 
degeneracy of the low laying energy bands is lifted, (a) One 
band exhibits four touching points between the valence and 
the conduction bands. Inset: points in /c-space where the 
valence and the conduction bands touch, (b) The other band 
becomes a Dirac cone. In these figures, the layer is unbiased 
and there is no ISO interaction. Other parameters have the 
values t = 1, t±_ = 0.2, a = 1, and k y — 27r/(v / 3«). 




►tR 

0.2 0.4 



FIG. 6: (Color online) (a) Zoom in on the K point of Fig.J^cQ 
and an intersection along the line k y = 27r/(VSa) of Fig. pi 
(b) Slope (a) of the Dirac cone as a function of the Rashba 
coupling tR. The scale of the a- axis depends on all parameters 
of the theory. Here: t = 1, t± = 0.2, a = 1, V = 0, A S o 0, 
and k y — 2tt/(v / 3g0- 



III. BILAYER GRAPHENE MODEL 



parameters, A so and £r, in units of t. Without SO inter- 
actions, we find the well known graphene spectrum with 
Dirac cones centered at the K and K' points in the recip- 
rocal spacef^see Fig.[2|a). The ISO interaction opens a 
gap, but respects the spin degeneracy of the energy bands 
(Fig.|2jb)). The Rashba term does not open a gap on its 
own, but it does lift the spin degeneracy, except at the 
k x = point, as it can be seen in Fig. [2|c). By zooming 
in on the region around the K (K') point, we see that the 
Rashba term splits the Dirac cones into four, as noted in 
Ref.|3l This behavior is depicted in Fig. [3] Note that this 
effect is missed in the approximation made by Kane and 
Mele, 2 which is effectively a zeroth order approximation 
of the Rashba term. If one takes linear terms in k into 
account, this effect is already present. Here, however, 
we keep the full expression for the spectrum, without re- 
sorting to approximations. The combined effect of the 
ISO interaction and the Rashba term breaks the particle 
hole symmetry (Fig. J2^d)). If the Rashba term is small 
(tR < A so), the gap is finite, 2 otherwise the gap closes 
(not shown). 



Before including the SO interactions into bilayer 
graphene, let us first consider the non-interacting Hamil- 
tonian. The form of the Hamiltonian of a multilayer 
graphenesystem depends on the stacking of the layers.^ 
For bilayers, however, there are only two possibilities. 
The lattice sites can lay exactly on top of each other, or 
they can be arranged in a Bernal stacking,^ in which the 
A sites of the upper layer (A\) lay on top of the B sites 
of the lower one (#2)? while the other sites (B\ and A2) 
lay opposite to a honeycomb center. We will assume the 
Bernal stacking here, because it is the most common one. 

The consequence of a Bernal stacking is that, in a first 
approximation, the only interlayer hopping is between 
A\ and B2 sites. It is straightforward to generalize the 
noninter acting monolayer Hamiltonian to a bilayer one. 
First, we define 

where the layers 1 and 2 are represented by the corre- 
sponding index. We introduce an interlayer hopping pa- 
rameter, t± ~ (0.1 — 0.2) t, and bias the bilayer system 
with a gate voltage V. This gate voltage can be tuned 
externally and is such that the lower layer has an elec- 
tric potential — V, while the upper layer has V. With 
this new parameters, the non-interacting Hamiltonian is 
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FIG. 7: (Color online) Energy spectrum of bilayer graphene for different values of the SO interactions. In this figure, the layer 
is biased, V = 0.25. Other parameters have the values t = 1, t± = 0.2, a = 1, and k y = 2tt / (\^3a). 



given by 



IV. INTRALAYER SO INTERACTIONS IN 
BILAYER GRAPHENE 



J d 2 k^(k)M° x8 *(k), 
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If we consider V = for the moment, we see in 
Fig. |4ja) that the spectrum is different from the one 
for monolayer graphene. The dispersion at the K (K') 
points, where the valence and the conduction bands 
touch, is no longer linear, but parabolic. This means 
that in a low energy approximation the quasi particles be- 
come massive (m ~ 0.054ra e )™ However, these particles 
are chiral with respect to sublattice pseudospin and are 
therefore massive chiral fermionsp^ which are new type 
of quasi particles, characteristic for bilayer graphene. 

If there is a nonzero voltage difference, a gap will open 
in the energy spectrum. In fact, not only a gap opens, but 
the parabolic bands are deformed i nto Me xican hats (see 
Ref. fT8l and the discussion in Sec. IV B). However, the 



energy bands remain spin degenerate, because a nonzero 
voltage on its own cannot lift the spin degeneracy; as we 
will discuss in the next section, interactions are required 
to reach this aim. 



A. No bias voltage 

Now that we have understood the single particle spec- 
trum for the bilayer graphene system, let us add SO in- 
teractions. We already saw that we have to distinguish 
between the ISO and the Rashba interactions, but for 
the bilayer system there is another subdivision, namely 
into intralayer and inter layer interactions. In this sec- 
tion, we will analyze the effects that the intralayer ISO 
and Rashba couplings have on the energy spectrum of 
the bilayer system. The effect of these interactions in the 
presence of a bias voltage are subsequently discussed. 

The ISO interaction respects the symmetries of a single 
graphene sheet. Since a graphene bilayer has a smaller 
symmetry group than a single layer, we expect this inter- 
action to be present in the planes of the bilayer system 
as well. The ISO interaction Hamiltonian is then given 

by 

H so = A so J d 2 k^(k)Mi? 8 <£(k), (7) 

(8) 
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Regarding the Rashba term, we expect that a perpen- 
dicular electric field gives rise to intralayer interactions 
in the same way that it did for a single sheet of graphene. 
Effectively, we have two copies of the monolayer Rashba 
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interaction, 



Hr = tR 



d 2 k^(k)M^(k), 



Mi 



M& 4 



Mg 4 



(9) 



(10) 



In Fig. [4j the energy spectrum of an unbiased bilayer 
graphene is shown for different values of A so and tR. For 
zero SO interactions, we observe the well known parabolic 
bands, which are spin degenerate (see Fig. Qa)). These 
degeneracies cannot be lifted by the ISO interaction on its 
own, which simply opens a gap in the spectrum (Fig. 4^b) 
and (c)). This is the same behavior as for monolayer 
graphene. Things become interesting when we consider 
the case of zero ISO coupling and a nonzero Rashba in- 
teraction (Fig.Qd)). The spin degeneracy of the bands is 
then lifted, but in a very particular way. In Fig. |5ja), we 
see that the K point splits into four points again, as for 
monolayer graphene. One of them remains at the former 
K point position and the others form a triangle around 
it, see the inset in Fig. |5|a). We will refer to these four 
points as the split K point. However, here very special 
features appear: Besides the two energy bands (conduc- 
tion and valence band) that touch at four points, there 
are two more bands that touch and form a Dirac cone, as 
it can be seen in Fig. (5^b). The center of this Dirac cone 
is exactly at the point in k space where, without Rashba 
interaction, the K point was located. This is also the lo- 
cation of the central of the four points that form the split 
K point. This central point is different from the other 
three. If we analyze Fig. |6|a), which shows a zoom-in 

of Fig. |4^d) and a cross-section cut of Figs. [5ja) and 
(b), we see that the central point, located at k x = 
27r/(3a), k y = 27r/(v / 3&) has, in addition to the Dirac 
cone, a linear crossing at very low energy, whereas the 
off-center points have only a higher order crossing. This 
is different from the case of monolayer graphene, where 
the K point splits into four equivalent points. The most 
striking feature of the intralayer Rashba coupling is the 
formation of a Dirac cone out of a parabolic band. If we 
perform a low energy approximation and use k = K + q, 
this Dirac cone has the dispersion relation E(q) = a\q\. 
The slope a, which corresponds to the velocity of the 
low energy excitations, depends on the parameters of the 



theory and cannot be determined analytically. We have 
plotted the slope of this cone as a function of tR for cer- 
tain parameter values in Fig. [6jb). 

If we set both tR and Aso unequal to zero (Figs. 4^e) 
and (f)), we see that depending on their relative values, 
a gap can open. For small ISO interactions, the gap 
stays closed, but the bands are heavily deformed in an 
asymmetric manner. Moreover, the particle-hole symme- 
try is lost. The split K point becomes so deformed that 
we cannot identify the four points any longer. If Aso 
becomes large enough a gap opens, but the asymmetry 
remains. 



B. Effect of a bias voltage 

If we add a bias voltage to a bilayer graphene system 
without SO interactions, the system becomes a semicon- 
ductor with a tunable gap. We will see in the following 
that SO interactions can heavily deform the energy dis- 
persion. 

As found earlier QSIMl we i3 Serve that for Aso = tR = 
the effect of the bias is to open a gap and to introduce a 
Mexican-hat-like shape in the lowest energy band around 
the K and K' points. The spectrum remains particle- 
hole symmetric and the bands remain spin degenerate, 
as is shown in Fig. [7^a). By keeping tR = and intro- 
ducing a finite Aso, the spin degeneracy is lifted, but 
only in the region around the K (K f ) points, as can be 
seen in Fig. |8ja). When Aso is not too large, the Mex- 
ican hat feature remains, but it disappears as Aso is 
increased further (see Figs. (7^b) and (c)). The gap at 
the K point (this is not the actual gap if V > Aso, be- 
cause of the Mexican hat structure), equals 2\V — Aso[, 
Therefore, during the transition from the Mexican hat to 
a parabolic band the gap closes at V = Aso, after which 
it opens again. Note that this behavior is similar to that 
described in Ref. 1, where both the ISO coupling Aso 
and a staggered sublattice potential X v can open a gap 
in monolayer graphene. If A so exceeds A v , a transition 
occurs between a normal insulator phase and a quantum 
spin Hall phase. Nonetheless, we have to be careful with 
this comparison, since a staggered sublattice potential 
in monolayer graphene is fundamentally different from a 
bias potential in a bilayer. For example, in monolayer 
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FIG. 9: (Color online) (a) Energy spectrum of a biased system, V — 0.25, with zero intralayer SO interactions and = 0.4. 
Compare with Fig. [TJa). (b) Energy spectrum of an unbiased bilayer system, with zero intralayer Rashba coupling, but with 
Aso = 0.15 and — 0.4. Compare with Fig. JZ^b). (c) Energy spectrum of an unbiased bilayer system, with zero intralayer 
Rashba coupling, but with Aso = 0.3 and t^ = 0.4. Other parameters are t=l,t± = 0.2, a = 1, and k y 2n/(V3a). 



graphene we do not observe a Mexican hat structure and 
the edge states are responsible for the phase transition. 
We do not take those into account here. In addition, 
although the ISO interaction alone cannot lift the spin 
degeneracy, it does so in the presence of a staggered sub- 
lattice potential in monolayer graphene. The lifting oc- 
curs for all values of k, except for k = 0, contrarily to the 
biased bilayer, where it is only significant around the K 
(K') points. Moreover, the effect is stronger in monolayer 
graphene than it is in bilayer graphene. 

As long as the Rashba coupling is zero, an analytical 
solution for the energy bands can be found. This solution 
is given by 



E(k) 



±,±,± _ 



± 
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72 



t 2 ± +2V' 2 +2t 2 \j\' 2 + 2A 2 so i] 2 
± (4+4t 2 ti| 7 | 2 + 16*V 2 | 7 | 2 
+ 16A| ^V ± 8t 2 ± A so Vv 



■ (11) 



It is clear that, if either V, Aso, or t± is zero the bands 
become degenerate, since in this case the last term 

with the ± sign vanishes. Another interesting fea- 
ture that we see only in a particular situation, where 
V, Aso 7^ and tR = 0, is a band crossing at k x = 
(see Figs. [TJb) and (c) and notice the inverted colors at 
the K and K' points), signalling that the k x — > —k x 
symmetry is lost. This band crossing can be seen an- 
alytically from Eq. (11). Note that |7/c| 2 is symmetric, 
while r] is antisymmetric under k x — > —k x . We conclude 
that, due to the linear term in 77, the energy bands satis- 
fies £;(fe) ± ' ± ' ± = Ei-k)^^. Therefore, the individual 
bands are no longer symmetric under this transforma- 
tion. 

Now, let us investigate the behavior of the system at 
finite £r, for Aso = 0- First of all, the spin degeneracy is 
lifted, except at k x =0. Second, the Mexican-hat feature 
evolves into something that looks more like an asymmet- 
ric farmers hat, see Fig.fJJcl) and Fig.[8fb). Note that if 
V and/or tR are increased, the asymmetry becomes more 
accentuated and the spectrum does not even look like a 



farmers hat anymore. We must emphasize that here we 
do not use any approximation for the energy spectrum, 
but we keep the full expression. This feature cannot be 
captured in a zeroth order approximation for the Rashba 
term, as used by Kane and Mele^, because at this order 
of the approximation the spectrum is symmetric around 
the K (K f ) point. 

If both Aso an d tR are finite, the spectrum becomes 
very complicated (Figs.JTJe) and (f)). As a general trend, 
Aso washes out the Mexican hat feature and at first, in- 
creases the difference between spin up and spin down 
bands around the K (K f ) points, although the spin de- 
generacy had been already lifted for all values, except 
at k x = 0, by the finite Rashba coupling. In addition, 
particle-hole symmetry is lost. There is no longer a band 
crossing, but the gap closes and opens again upon in- 
creasing Aso- Depending on the parameters, a situation 
can occur where the bands do not touch, but they have 
common energies, thus there is no gap in the system (see 
Fig.§c)). 



V. INTRA- AND INTER-LAYER SO 
INTERACTIONS IN BILAYER GRAPHENE 



In the previous section we accounted for tunneling be- 
tween the two layers, but considered only intralayer SO 
interactions. Now, we investigate the effect of SO inter- 
actions between lattice sites in different layers. Since 
the ISO interaction depends on the symmetry of the 
graphene plane, it is not obvious if there is any inter- 
plane ISO interaction at all. Hence, we focus on the 
Rashba term. This term only exists in graphene mono- 
layers if the z — > — z mirror symmetry is broken, for ex- 
ample by a perpendicular electric field. Furthermore, the 
Rashba coefficient can be tuned by varying this electric 
field. Therefore, we consider here a bilayer system in the 
presence of a tilted electric field. The inplane component 
of the electric field (En) gives rise to an interlayer Rashba 



coupling that is a generalization of Eq. Q, 

Hk = -itk E a h ( s x s ) • E|| &i,2 + h.c, (12) 

E|| = (cos(/>, sin^, 0) T , 

where we have chosen to absorb the magnitude of the 
electric field already in the constant t^. The orientation 
of E|| is determined by 0, but the results will be inde- 
pendent of (j) and therefore we choose = arbitrarily. 
The unit vector connecting the two lattice sites and 
Bi 2 is given by — z and this explains the minus sign in 
comparison with Eq. [4j In /c-space, Eq. ( 12 ) becomes 
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The effect of this interlayer Rashba interaction depends 
heavily on the other parameters in the theory. Without 
any intralayer SO interactions, the result of a nonzero 
tjz is effectively a modification of the interlayer hopping 
parameter, 



t±-*xVi + (*fl)7*i- 



(13) 



The effect is a slight deformation of the energy bands. 
Only if tjz > 0.3 the shift becomes significant. For a 
biased system, tk will flatten the Mexican hat, as it can 
be seen in Fig. [9^a). 

The effect of interlayer Rashba coupling is most visi- 
ble in bilayer systems with zero bias, but with intralayer 
ISO interactions. Although the energy spectrum can 
be solved analytically, the equations become too com- 
plicated to handle. However, it is clear that the effect 
is more than a shift of the interlayer hopping parameter. 
The spin degeneracy of the energy bands is lifted and we 
see a Mexican hat feature appear in the low laying energy 
band (see Fig.^b)). These low energy bands are shifted 
towards the Fermi level and as a consequence the gap 
between the valence and the conduction bands becomes 
smaller. If A so is small, the spin degeneracy of the low- 
est energy band stays intact at the K point, but is lifted 
around it. However, if Aso increases, this degeneracy is 
lifted and eventually shifted to a degeneracy between the 
two conduction bands that lay in the middle, as it can be 
seen in Fig.|9jc). In comparison with the bilayer without 
intralayer interactions, the effect of nonzero manifests 




FIG. 10: (Color online) (a) Energy spectrum of an unbiased 
system with zero intralayer ISO interaction, but tR = 0.2 and 
= 0.4. Compare with Fig. |4[d) . (b) Zoom in on the K 
point of (a). Compare with Fig. |6[a). Other parameters are 
t=l,t±= 0.2, a = 1 and k y = 27r/( v / 3a). 



itself also for small values of this parameter. In a bi- 
ased system with intralayer ISO interactions, the effect 
of the interlayer Rashba coupling is visible, but its influ- 
ence becomes less important as the bias becomes larger. 
The effect of a nonzero in this case is to increase the 
splitting of the bands around the K (K 1 ) points and if 
V > Aso, the Mexican hat is flattened (not shown). 

Let us now consider an unbiased system with zero 
ISO coupling, but with nonzero interlayer and intralayer 
Rashba coupling. We would have such a system if a tilted 
electric field is present. The relative strength of both in- 
teractions can be tuned independently by changing the 
parallel and perpendicular components of the electric 
field. This case is a generalization of Fig.EFd), which 
was described in more detail in Figs. [5] and [6^ The most 
striking effect of a nonzero is the lifting of the spin 
degeneracies at the K (K f ) points, with as main conse- 
quence the destruction of the Dirac cones (see Fig. 10). 



The way the lowest laying energy band splits is also ef- 
fected by the interlayer Rashba coupling. Recall that the 
K (K 1 ) point splits into four due to intralayer Rashba 
coupling and that the central point had a linear crossing 
at very low energy scales. This linear crossing is modi- 
fied to a higher order crossing by the interlayer Rashba 
interaction, as it can be seen in Fig. 10 Note that the 
energy scale of the zoom in at the K point in Fig. [To|has 
decreased by a factor fifty with respect to Fig. [6] This is 
also due to the interlayer Rashba coupling. 

If all parameters are nonzero the spectrum is very com- 
plicated. Interlayer Rashba coupling does deform this 
spectrum, but we could not detect any special feature 
that would justify exhibiting them. 



VI. CONCLUSIONS 

We studied a graphene bilayer including both, the in- 
trinsic and the Rashba SO interactions and we found that 
these interactions can modify the energy dispersion in a 
non trivial way. 

First, we concentrated on the unbiased system and con- 
sidered only intralayer interactions. We observed that 
ISO interactions can open a gap in the system, in the 
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same way as they do for monolayer grapheneP On the 
other hand, if only an intralayer Rashba SO interaction 
is present, the energy bands are completely different than 
for monolayer graphene. The K (K f ) points still split into 
four, but at very low energies the off center points have 
now a higher order crossing instead of the linear cross- 
ing present in monolayer graphene (Fig.^a)). The most 
striking feature, however, is the formation of a Dirac cone 
out of an energy band that once was parabolic (Fig.[5jb)). 
This Dirac cone is located exactly at the former K (K') 
point. The intralayer Rashba interaction not only lifts 
the spin degeneracy of the energy bands, but also changes 
their individual behavior. If both the intralayer ISO and 
Rashba terms are nonzero, the Dirac cone is destroyed 
and the particle- hole symmetry is broken. 

The presence of a Dirac cone in unbiased bilayer 
graphene is a very special feature. Together with the split 
K (K 1 ) point, the spectrum that we found will give rise 
to two different low energy excitations, one of which is 
massless. The speed with which this massless excitation 
travels depends on the Rashba constant (see Fig. [6jb)), 
which can be tuned by the perpendicular electric field. 
We expect that the dispersion relation shown in Fig.|6ja) 
could be observed with ARPES experiments, which have 
already successfully demonstrated the Dirac dispersion 
in monolayer graphene!^ 

In a biased system with nonzero ISO interactions, the 
spin degeneracy of the energy bands is lifted. The split- 
ting occurs only around the K (K f ) points. We also 
observe a band crossing at the k x = point (Fig. [7jc) 
and (d)). A bias voltage, in combination with the in- 
tralayer Rashba coupling, destroys the Dirac cones and 
the spectrum becomes asymmetric around the K (K') 
points. 

Next, we considered an interlayer Rashba interaction 
between the planes. This interlayer interaction would in 
principle be expected to be small compared to the in- 
tralayer one because of the larger interlayer atomic sep- 
aration. However, this effect could still be important if 
pressure is applied to approach the two layers. In our 
model, the interlayer Rashba coupling finds its origin in 
the presence of a tilted electric field. In a bilayer sys- 
tem with no intralayer SO interactions, this interaction 
causes effectively a shift of the interlayer hopping pa- 
rameter. However, in a system where the intralayer SO 
couplings are nonzero, we see a clear effect in the energy 
spectrum. In an unbiased system with nonzero Ago, the 
spin degeneracy of the bands is lifted around the K (K f ) 
points and a Mexican hat feature appears. The Mexi- 
can hat feature is known to arise in the bilayer system in 
the presence of a bias volt age Here, however, we found 
that it can also appear without a bias, but solely due to 
SO interactions. Moreover, the Mexican hat in the en- 
ergy dispersion is fully spin polarized. If A so becomes 
large enough, the degeneracy of the two lowest laying 
conduction bands at the K (K 1 ) points is shifted to the 
two middle bands. The system then becomes isospin de- 
generate. Indeed, in the presence of tunneling, a bilayer 



can be described as a two-level system, where the energy 
bands in each layer have combined into symmetric and 
antisymmetric energy bands, separated by a gap given 
by the tunneling energy. If we represent the asymmet- 
ric band by an isospin up and the symmetric one by an 
isospin down, we see that the ISO interaction can lead to 
an isospin degeneracy at the K and K' points, although 
the spins remain fully polarized. 

If the ISO interaction is absent, but the intralayer 



Rashba is nonzero, we have seen in Fig. 10 that the ef- 
fect of the interlayer Rashba interaction is to destroy the 
Dirac cone at the K (K f ) point and to modify the way the 
K (K r ) point splits into four. We can no longer observe a 
linear crossing for the central point and the energy scales 
associated with this splitting are substantially smaller. 

Finally, we would like to discuss the possibility of ob- 
serving experimentally the effects that we have described 
above. In monolayer graphene, the current estimates 
are that the ISO interaction is very small (0.0011 — 0.05 
meV). 6 However, it is already possible to tune the Rashba 
coupling in a graphene layer on a Ni substrate up to 
tR ~ 0.2 eV.^ Because these values should be represen- 
tative for bilayer graphene as well, we expect that the re- 
sults we found involving the Rashba interaction are well 
within the experimental reach. We think it should be 
possible to detect the Dirac cone that arises from the in- 
tralayer Rashba term in an unbiased bilayer. If this Dirac 
cone can be detected and if it is destroyed by an in-plane 
electric field, we would have an indication that indeed an 
interlayer Rashba interaction is present in the system. 

The values we used for Aso are larger than indicated 
experimentally. However, we should recall that the same 
kind of system can be engineered using cold atoms in 
optical lattices, and in this case there is much less con- 
straint on the parameters of the model. Our aim here was 
mainly to draw a comparison of the different effects of 
the Rashba and ISO interactions to determine the trend 
introduced by each one. 

A next step would be to include edge states in the 
model. If we regard the intralayer ISO interaction and 
the bias voltage in bilayer graphene as being compara- 
ble with the ISO interaction and a staggered sublattice 
potential in monolayer graphene, there is a possibility 
that the bilayer system would exhibit a phase transition, 
equivalent to the one described in Ref [U It is already 
known that there are two types of edge states in bilayer 
graphenep^ but SO interactions have not yet been taken 
into account. We hope that our results will motivate 
further theoretical studies and experiments in the field. 
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